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Abstract. In this paper we consider the field equations for linearized gravity 
and other integer spin fields on the Kerr spacetime, and more generally on 
spacetimes of Petrov type D. We give a derivation, using the GHP formalism, 
of decoupled field equations for the linearized Weyl scalars for all spin weights 
and identify the gauge source functions occuring in these. For the spin weight 
Weyl scalar, imposing a generalized harmonic coordinate gauge yields a gen- 
eralization of the Regge- Wheeler equation. Specializing to the Schwarzschild 
case, we derive the gauge invariant Regge- Wheeler and Zerilli equation directly 
from the equation for the spin scalar. 
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l/-) ■ 1. Introduction 

GN ■ In his 1965 paper [37] Penrose showed that, at least formally, all solutions of the 



massless spin s field equation on Minkowski space can be obtained from solutions 
of the spin (s — |) equations. By repeating this process 2s times, one finds that 
solutions of to the spin equation (i.e. the free scalar wave equation) are potentials 
for the massless spin s field. This is an example of the spin-raising and -lowering 
, transformations discussed in more detail in [?TJ section 6.4], see also [411 section 

6.7]. 

The following special case of the construction is relevant here. Given a 2 index 
Killing spinor, i.e. a spinor Kab solving the equation 

V ( /'Ab C) = 0, 

and a solution 4>abc...d = 4>{ABC...d) °f f ne s P m s zer0 rest- mass equation 

V AA '4>ABC...D = 0, 

the spin-lowered field (j>c...D = 4>abc...dK ab is a solution of the spin (s — 1) zero 
rest-mass equation. 

The analysis of linear field equations on Minkowski space is a key step in the proof 
of the non-linear stability of Minkowski space [111 IT2"] . For the case of Minkowski 
space, the linearized Bianchi equation is precisely the massless spin 2 equation. 
Thus, the discussion of Penrose shows that the linearized stability problem for the 
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Einstein equation on Minkowski space can essentially be reduced to a study of the 
scalar wave equation. 

The problem of non-linear stability of Minkowski space, solved in [TJ], see also 
|32| . can be viewed as a warm-up for the black hole stability problem, i.e. the 
problem of proving the non-linear stability of the Kerr spacetime in the class of 
asymptotically flat vacuum spacetimes, see e.g. [UE3H2], an< ^ references therein. 
The black hole stability problem adds several levels of difficulty over the problem of 
stability of Minkowski space. The aspect which we shall focus on in the present dis- 
cussion is that the background spacetime (or more properly stated, the asymptotic 
state for the evolution) fails to be conformally flat and consequently the relation 
between the spin wave equation, viewed as a model problem for the full non-linear 
stability problem, and the equations of linearized gravity fails to be as close as on 
the Minkowski background. 

In particular, the equation for linearized gravity (i.e. the linearized Bianchi 
system) on a non-conformally flat vacuum background is not the spin-2 system, 
but has a non-trivial right hand side. In fact, the massless spin 2 equation in a 
spacetime of Petrov type D has only trivial solutions, cf. [5]. It follows that the 
spin-lowering and -raising transformations cannot be applied directly to Maxwell 
or linearized gravity on a background spacetime which is not conformally flat. 

However, if we consider the field equations of integer spin on a vacuum type D 
spacetime, of which Kerr is a special case, then it turns out that an analogue of 
the spin-lowering transformation does yield useful results, even for the equations 
of linearized gravity. This idea was alluded to already in the paper of Jeffryes, see 
[2"91 p. 340]. Although not stated as explicitly, related ideas play an important role 
in the work of Fackerell and Crossman [Ml [18] . 

From the point of view of the black hole stability problem, the two-parameter 
Kerr family of rotating black hole solutions are the vacuum type D spacetimes of 
most immediate interest. However, most of the results in this paper are valid for 
general vacuum type D spacetimes. These include the Kerr-Taub-NUT spacetimes, 
see [33 section 21.1], see also |16| . 

1.1. Spin- lowering the Maxwell field. RecalQ [5U] that any vacuum type D 
spacetime admits a Killing spinor of the form 

K AB = y- 1/3 (A L B y (1.1) 

A Killing spinor which in addition satisfies the condition 

S7 A , B K AB +\7 A B 'K A , B , =0 

corresponds via 

K a b = i{K AB e A > B > — K A , B >e AB ) 

to a Killing- Yano tensor, i.e. a skew 2-tensor K a b = Kwm satisfying the Killing- 
Yano equation 

V ia K b)c = 0, 

see [351 031 130] for further information. 

Let 4> AB be the Maxwell spinor, i.e. a solution of the massless spin 1 equation. 
Then (f> AB K AB = ^ 2 1/3 ^, where fa is the spin weight Maxwell scalar and this 
rescaling of <pi solves a wave equation with potential 

(□ + 2* 2 )(*2 1/3 0i) = O, (1.2) 



Here and below we use the conventions and notations of the Newman-Penrose (NP) and 
Geroch-Held-Penrose (GHP) null-tetrad based formalisms, and the associated two-spinor formal- 
ism, see [36] [25] [40] [41] and section [2] below. 
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where □ = V a V a is the d'Alembertian of the spacetime and ^> 2 is the spin weight 
Weyl scalar of the background. This is precisely the wave equation derived for 
the rescaled spin- weight Maxwell scalar on Kerr by Fackerell and Ipser [19] . who 
also argued that 4>i can be used as a potential for the Maxwell field on the Kerr 
spacetime. 

The special case of this spin weight wave equation for the Maxwell field on 
the Schwarzschild spacetime was used recently in the proof by Blue [5] of decay 
estimates for the Maxwell equation. In this work, Blue was inspired by Price 02] 
who showed that the Regge- Wheeler [43] wave equation for axial perturbations of 
the Schwarzschild spacetime can be viewed as a wave equation, with potential, for 
a rescaled version of the imaginary part of the spin weight linearized Weyl scalar. 

1.2. Gauge invariant equations for linearized gravity. We follow the conven- 
tion of Price and others in discussing perturbations of GHP quantities and let a 
subindex A denote quantities defined on the background and indicate first order 
perturbed quantities with a subindex B. See |47l 17] for detailed treatments of per- 
turbation theory in the context of tetrad based formalisms. The Regge- Wheeler 
wave equation in the form derived by Price can be written as 

(□ + 8^ 2A )(* 2 ^ /3 3m* 2S ) =0, 

on the Schwarzschild spacetime, where in Schwarzschild coordinates ^>2A = —Mr~ 3 . 

The approach of Price does not generalize directly to include polar perturbations. 
These had previously been treated by Zerilli [52], see also Moncrief [35] . who de- 
rived a non-local wave equation governing a gauge invariant potential for the polar 
degrees of freedom, and related work by Bicak [3] on perturbations of the Reissner- 
Nordstrom solution. Among the difficulties in generalizing Price's approach to the 
Regge- Wheeler equation to cover general perturbations even of Schwarzschild is 
that in the background the ^2A is non-zero and hence ^2B fails to be gauge invari- 
ant. In fact, it is only the linearized Weyl scalars ^ob, ^4B of extreme spin weights 
2, —2, which are both coordinate and tetrad-gauge invariant, and hence it is only 
those which can be directly viewed as physically measurable quantities. 

Teukolsky [151 HH] showed, by calculating in the NP formalism on a Kerr back- 
ground, working in the principal Kinnersley tetrad, that suitably rescaled versions of 
the extreme spin weight Maxwell (s = 1, —1) and linearized Weyl scalars (s = 2, —2) 
on the Kerr spacetime satisfy decoupled and separable wave equations. The result- 
ing system is usually called the Teukolsky Master Equation (TME). 

Shortly after the work of Teukolsky, Ryan [U] showed that the vacuum Teukolsky 
system for linearized gravity can be derived simply by projecting the Penrose wave 
equation, i.e. the covariant tensor wave equation 

ORabcd = RabefRcd Cf + 2(R a ec f R^ d ~ Raedf R^ c ) , 

satisfied by the curvature tensor, on a principal null tetrad, and linearizing. This 
theme has been taken up and generalized to arbitrary vacuum backgrounds by Bini 
et al. [HH], where all equations for the ^'s and 4>'s are calculated from components 
of a generalized de Rham operator acting on the Riemann and Maxwell tensor. 

Now consider a general vacuum type D background spacetime. Working in a 
principal tetrad, let ip s be one of the fields * b, * 2 ^ /3 *4b,0 o , * 2 a /3 02 (letting s 
be the spin weight of the field). Further, let 

B a = -(pn a - Tfh a ). 

Define a generalized wave operator acting on properly weighted quantities by 

B p = g ab {Q a + P B a ){Q b + P B b ), 
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where O a is the weighted GHP covariant derivative, see section [2] for details. Then, 
cf. [1], the TME S for spin s fields can be written in the form 

(Q 2s - AsH 2A )^ s = 0. (1.3) 

As shown by Teukolsky, for the Kerr case, equation (|1.3p can be separated into 
radial and angular equations. Due to this fact, it has been possible to formally 
analyze the Teukolsky system and its solutions and many interesting discoveries 
have been made. Among these are the relation of the separability of the system 
to the presence of the Carter constant and corresponding symmetry operators, the 
calculation of the separation constants, as well as the proof of mode stability for 
the Teukolsky system |51) . 

However, it is not clear from these works that the Teukolsky system is well suited 
for the analysis of the asymptotic decay properties of the higher spin fields. Among 
the difficulties encountered in attempting to analyze the Teukolsky system are the 
facts that it has a long-range potential and lower-order terms with slowly decaying 
coefficients. See [S5| and citations therein for discussion. 

1.3. Spin lowering the linearized Weyl field. Recall that on a vacuum type D 
spacetime, the spin-lowered Weyl field i^abcdK , where Kab is a Killing spinor, 
satisfies the Maxwell equation, see 1311 §3.8], see also [33] • The same statement 
holds for the linearized Weyl spinoio {S^abcd on Minkowski space, and further 
in that case lowering the spin by 2 gives {8iP)abcdK ab K CD which is a solution to 
the spin wave equation. We now consider the equations satisfied by these fields 
on a vacuum type D background. 

Let (Sip) abc d be the linearized Weyl spinor on a vacuum type D background, 
and let Kab be the Killing spinor as in Then, expanding the spin 1 field 

(SiP)abcdK cd which corresponds to a skew 2-tensor, into weighted scalars gives 

— 1/3 

the rescaled linearized Weyl scalars 4>t-\ — ^2A ^iB, i = 1,2,3, of spin weights 
1, 0, —1 and additional terms arising from linearized tetrad. Thus, by analogy with 

~ 1/3 2/3 ? 

the above, it is reasonable to suppose that the fields <fto = ^ 2 A ^2A = 

^2A^3B of spin weights 1,-1, respectively, satisfy an analogue of the Teukolsky 

* 2 /3 

system TME S , for ,s = 1,-1, while the spin weight field (f>i = ^ 2 A *2B can be 
expected to satisfy an analogue of the Fackerell-Ipser equation. Further, since the 
non-extreme linearized Weyl scalars fail to be gauge invariant, one also expects the 
corresponding equations to contain gauge potential terms. 

For the spin weight one case, a calculation, see section l3~2l shows 

(S 2 - 4* 2 )(*2i /3 *iB) = -6*2A t0>' + V - ~P>b - (8' + 2r' - f)a B + 2*ib], 

The equation satisfied by *f> 2A&3B is similar. The right hand side of the equation for 
VPib corresponds to the right hand sides of the equations for 0qb on a charged type 
D background with ^2A playing the role of the spin weight Maxwell scalar 4>ia, 
see section l3~2l The analogous statement holds for after applying a prime. As 
we shall see, the conditions that the right hand sides of the equations for ^ib, 
are zero are tetrad gauge conditions. For the case of the perturbed Maxwell field 
this corresponds to turning off the background charge, i.e. imposing the condition 
4>iA — 0. The fact that the just mentioned tetrad gauge conditions can be viewed 
as the "ghost" of the background charge motivated Chandrasekhar [TOj p. 240] to 
use the term phantom gauge. In fact, Chandrasekhar showed that by imposing the 
phantom gauge, the rescalings ^> 2 ^ 3 ^ib, ^2~a^3B satisfy the TME S , for s = 1, — 1, 
respectively. 



Here we use a <5 to denote linearized quantities, eg. (S^P)abcd in order to avoid confusion 
with spinor indices. 
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We show in section|3]below that the phantom gauge can be viewed as prescribing 
a gauge source function for the tetrad degrees of freedom along the lines of Friedrich 
|23| . with the linearized Weyl field itself as part of the gauge source. The phantom 
gauge was studied by Chandrasekhar from a formal point of view only, and the 
possible implications of this procedure for the hyperbolicity and well-posedness of 
the linearized Einstein equations were not analyzed in his work. We show here that 
the phantom gauge condition is compatible with a well posed system of equations 
for linearized gravity. 

Finally, we consider the spin weight linearized Weyl scalar. As the background 
^2A is non- vanishing, one has that ^2B is coordinate gauge dependent (but tetrad- 
gauge independent). Motivated by the previous discussion we consider the equation 
satisfied by the spin weight field ^ 2 A obtained by lowering the spin of the 

the linearized Weyl field by 2. A calculation, cf. section shows that this rescaled 
spin weight linearized Weyl scalar solves the equation 

(□ + 8* 2 )(* 2 j /3 * 2B ) = -3D B *^ 3 , (1.4) 

where the right hand side is the first order perturbation of the wave operator, 
acting on the background spin weight zero Weyl scalar. In section [5] we show that, 
restricting to the Schwarzschild case, equation (|1.4|) contains all of the information 
in the Regge- Wheeler and Zerilli-Moncrief systems by giving a direct derivation of 
these systems starting from (|1.4[) . 

The condition that the right hand side of (| 1 .4[) vanishes can be viewed as a gener- 
alized harmonic coordinate condition. It is worth noting that this gauge condition 
can be imposed also in the Schwarzschild case. Imposing this generalized harmonic 
gauge condition, the spin weight zero linearized Weyl scalar satisfies the scalar wave 
equation 

(□ + 8*2A)(*2"i /3 *2fl)=0, (1.5) 

which can be viewed as a generalization of Price's version of the Regge- Wheeler 
equation not only to the full set of perturbations of Schwarzschild but also to 
perturbations of Kerr. Lun and Fackerell |33| considered the situation on Schwarz- 
schild and argued formally that by imposing a suitable gauge condition, one obtains 
equation (|1.5p (specialized to the Schwarzschild case). 

We further point out that a generalized harmonic gauge condition with a gauge 
source function involving can also be used to modify the potential in (|1.5[) 
so that the equation becomes the Fackerell-Ipser equation. However, as already 
discussed by Crossman and Fackerell [HI [TS] this is possible only in the rotating 
case, and in particular for the Kerr family of spacetimes involves a division by a. 
As in the case of [33] , the discussion in the papers [T3J [TH] is quite formal and the 
gauge conditions are there not expressed in terms of gauge source functions. It is 
interesting to note that the just mentioned work of Crossman and Fackerell took 
as a starting point the Maxwell equation for the spin-lowered Weyl field in a type 
D spacetime, and its linearization. This has been carried further in the work of 
Ferrando et al. |20| where gauge conditions yielding an exact Maxwell system for 
the linearized, spin-lowered Weyl field have been considered. 

Now, following the approach taken by Blue for the case of Maxwell on Schwarz- 
schild, where the spin weight zero scalar was used as a potential for the Maxwell 
field, it is an interesting possibility to use ^ 2A *2B as a potential for the full 
linearized gravity system. This would allow one to reduce decay estimates for lin- 
earized gravity on Kerr to the study of the scalar wave equation with potential 

As is indicated by the discussion above, there is a great deal of freedom in using 
gauge conditions to change the nature of the (tetrad based) linearized Einstein 
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equations. It is to be expected that this remark applies equally to the full, non-linear 
system of Einstein equations. The implications of this remain to be considered. 
One could in principle go further, and make use of the gauge dependence of ^2B 
to remove the potential from the equation and achieve a setup (in the rotating 

_ 2 / 3 2/3 

case) where *& 2 a ^ 2B s °l ves the scalar wave equation □(* 2y J /J * 2 s) = 0. We 
point out that the gauge conditions chosen by Chandrasekhar, cf. [101 §82], in his 
considerations of linearized gravity on the Kerr background included the conditions 
^fiB = ^2B = ^35 = 0. This type of gauge condition will not be considered in 
detail here. 

For linearized gravity on Kerr, among the questions which should be considered 
are the choice of potential for linearized gravity and the field equation governing 
this. For the gauge invariant scalars ^ob,^ab these issues have been extensively 
discussed in the literature. Making use of gauge conditions as discussed above 
opens up interesting new possibilities. However to make full use of these, the 
problem of reconstructing the full solution of linearized gravity from e.g. ^2B must 
be considered. 

1.4. Overview of this paper. The plan of this paper is as follows. In section[2]we 
set up notation and give a brief overview of the GHP formalism and its specialization 
to vacuum type D backgrounds. We also discuss there gauge issues that arise when 
working in a tetrad based formalism. Section [3] starts by introducing the properly 
weighted generalized wave operators which occur in the GHP formalism and in 
the Teukolsky system. Further, we give there a derivation in the GHP formalism 
of the equations for linearized gravity, cf. section 13.21 The gauge nature of the 
non-trivial right hand sides of the equations for the scalars of non-extreme spin 
weights is discussed in section 0] where we also give gauge-fixed versions of these 
systems which lead to new potentials for linearized gravity satisfying well-posed field 
equations. In section [5] the gauge invariant Regge- Wheeler and Zerilli equations for 
Schwarzschild background are derived from (|1.4D . 

2. Preliminaries and notation 

We use the conventions and notations of |25| . In particular we use abstract index 
notation with lower case latin indices for tensors and upper case latin indices for 
spinors. For tetrad indices we use lower case fraktur font, while for coordinate 
indices we use lower case greek letters. Unless otherwise stated we shall consider 
only vacuum spacetimes (M,g ab ) of dimension 4, with signature H . 

2.1. GHP formalism. The Geroch- Held- Penrose (GHP) null-tetrad formalism 
|25| allows one to represent the Einstein equations in a compact form, and gives an 
efficient tool for calculations. Since we will make heavy use of this formalism and 
its properties, we give, in order to make the paper reasonably self-contained, a brief 
description of its main features. 

Consider a null tetrad (e„) = (l a , n a , m a , rh a ) consisting of two real null vectors 
l a ,n a and two complex linear combinations of spatial vectors m a ,rh a , normalized 
such that the only non- vanishing inner products of the tetrad vectors are 

l a n a = —m a m a = 1. 

In order to avoid clutter we suppress the abstract index on the tetrad vectors e a 
when convenient. The coframe e° is defined by the relations e a (e^) = 5%. We note 
the useful relations 

9ab = l a n.b + nJb - m a fh b - rh a m b , (2.1a) 
5 a b = l a n b + n a l b - m a fh b - m a m b . (2.1b) 
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A choice of a null tetrad picks out a 2-dimensional subgroup of the Lorentz group in 
each tangent space (and hence also a reduction of the principal SO+(3, 1) bundle of 
(M, g), see below) which preserves the null planes spanned by l a , n a and the spatial 
planes spanned by m a , fh a . These can be represented in terms of a non- vanishing 
complex field A by the boost rotations 

l a -> AA>, n a -> A _1 A _1 n a , (2.2a) 

and the spin rotations 

m a -> AA _1 m a , m a -> A _1 Am a . (2.2b) 

Projecting tensor fields on the spacetime on the null tetrad gives a representation of 
these fields in terms of collections of tetrad components which are simply complex 
fields on the spacetime. In general, a scalar field r\ defined by projecting a tensor 
field will transform as 

r) -> \ p X q rj , 

for some integers p, q, under the above defined action of A. A quantity rj which 
transforms according to the above rule is said to have type {p, q} and fields with 
well defined type are referred to as weighted quantities. Note the notion of weighted 
quantity extends to tensors. In particular, the tetrad vectors l a , n a , m a , rh a have 
types {1,1}, {—1,-1}, {1,-1}, { — 1,1}, respectively. It is useful to note that 
the type is additive under multiplication and hence the weighted quantities form 
a graded algebra. The spin and boost weights s,r of a weighted quantity r\ are 

• s = |(p-?)) r = \{p + i)- 

The following formal operations take weighted quantities to weighted quantities, 

"(bar) : l a -> l a , n a -> n a , m a -> m°, fh a -> m a , {p,q} -> {g,p}, 

'(prime) : l a -> n a , n Q P, m a fh a , m a -> m a , {p, q} -> {-p, -q}, 

*(star) : l a -> ?n a , n a -> -m Q , m a -> -/ a , ?ri a -> n a , {p, g} -> {p,-g}. 

(2.3) 

The bar operation acting on a weighted quantity is simply the complex conjugation 
of the field. We have fj = n and 77" = (— l) p+q n (note however that we shall consider 
only fields with p + q even). Further, the bar ~ and prime ' operations commute, 
while the star * operation commutes with neither of these. Thus, the star operation 
has to be treated separately from the bar and prime operations. 

For the case of an orthonormalized tetrad, the Levi-Civita connection can be 
represented in terms of 24 independent real connection coefficients. In terms of 
the null tetrad introduced above, the connection coefficients T c ab = e^(V ea e^) com- 
bine into 12 complex scalars, called spin coefficients. Of these only 8 are properly 
weighted, and can be represented in terms of the quantities 

K = m h rV a lb, o- = m b m a VJb, P = m b fh a VJ b , r = m b n a VJ b , (2.4) 

and their primes a', p', t'. The effect of the star operation on the spin coefficients 
and their primes and complex conjugates can be calculated directly from (|2.4|) . or 
see [23 p. 878] for a list. The types of the spin coefficients are 

«:{3,1}, a: {3,-1}, p:{l,l}, r:{l,-l}, 

and the types of their primes are given according to (|2.3[) . 
The remaining 4 spin coefficients 

13 = ^{n b m a V a l b - m b m a V a m b ), 
e=^(n b l a V a l b -m b l a V a m b ), 
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and their primes /?', e', are not properly weighted and are not used explicitly in the 
GHP formalism. 

Let C* denote the non-zero complex numbers. As discussed by Ehlers |17) . 
see also [26l [27] , the choice of two null directions gives a reduction of the principal 
SO+(3, 1) frame bundle of the spacetime to a principal C» bundle B with the action 
of C* z : (A, e„) t— > e a .z(X) given by (|2.2p . The weighted quantities may be viewed 
as sections of associated complex line bundles! £ {p ' q} = B y. z { v .q} C determined by 
the representation z^ p ' q ^ : (A, v) — > X p X q v of C* on C. The restriction of the Levi- 
Civita connection to the reduced bundle B induces a connection on the bundles 
£{p,q} given by 

where w a is the connection form 

uj a = -('la + en a + f3'm a - /3m a 

= ^{n b V a h + m b V a fh b ). 

Under a gauge transformation u a transforms as 

, V a A 

UJ a .Z = LU a H — , 

Note that oj' a = —LJ a - It follows that 

for any properly weighted quantity rj. 

The GHP operators p,p',3, 5' are defined as the a covariant derivative along 
the tetrad vectors, 

b = re a , b' = » a e a , 8 = m a e a , 3' = m a e a . 

The action of the bar, prime and star operations on the GHP operators follows from 
their action on the tetrad vectors. Expanding 8 a in terms of the GHP operators 
gives 

©a = laV + n a t> -m a Q' — ffl a 8 . (2.6) 

In terms of the graded algebra of weighted quantities, the covariant derivative a 
as well as the GHP operators satisfy a (graded) Leibniz rule. Further, as remarked 
above, the notions of properly weighted quantity extend to differential forms and 
more general objects, and hence also the weighted covariant derivative lifts to act on 
such objects. In particular, the tetrad elements themselves are properly weighted 
quantities, and the action of O a and the GHP operators on these can be read off from 
the definitions. We have for example bZ c = 6 C l a Q a lb. Expanding this out using 
(|2.1b|) . (|2.5[) and the definitions of the spin-coefficients gives bZ c = ~Rm c — nfh c . 
The following equations 

bl a = -nm a — Kfh a , b m a — —f'l a - Kn a , (2.7a) 

b'L = "TO,, - ™ o , b' m a = — «'i a - ™a, ( 2 -7b) 

bla — —pm a — afh a , bm a = -a'l a - an a , ( 2 -7c) 

b'l a = -am a - prh a , 8'm Q = —p'l a - pn a - (2.7d) 

and their primes and complex conjugates give the complete set of relations. 

The Leibnitz rule for covariant derivative together with (|2.7|) allows one to read 
off the action of the GHP operators on tensors projected on any combination of 
tetrad vectors, and hence covariant tensor equations can be expressed equivalently 
as collections of scalar equations in the GHP formalism. 



By the same construction we can also treat e.g. differential forms as weighted quantities. 
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The 10 degrees of freedom of the Weyl tensor can be represented by the five 
weighted Weyl scalars 

*o - W abcd m a l b l c m d , *! = W abcd n a l b l c m d , * 2 = W abcd m a l b fh c n d , 

* 3 = W abcd n a l b rh c n d , <J 4 = W abcd rh a n b n c fh d . 

Similarly, the Maxwell field strength can be represented by the 3 Maxwell scalars 

fa = F ab l a m b , fa = ^(F ab l a n b + F ab fh a m b ), fa = F ab fh a n b . 

We shall refer to the spin coefficients and the Weyl and Maxwell scalars collectively 
as GHP quantities. 

The Weyl scalars i = 0, ...,4 have types {4 — 2i,0} while the Maxwell 
scalars fa, i = 0, 1, 2 have types {2 — 2i, 0}. The prime operation gives ^ = ^>4-i, 
i = 0,..., 4, and $ = -fa-i, i = 0, 1, 2. Further, <F* = ^ and fa* = fa. 

We now state, modulo prime and star operations, the Einstein, Bianchi and 
Maxwell equations in GHP notation, specialized to the vacuum case. Working in 
a tetrad formalism, the Einstein equation takes the form of a system of first-order 
equations for the connection coefficients, i.e. in the GHP setting, for the spin 
coefficients. These are given by 

8p - 8'tr = (p - p)T + (p' - P ')k - (2.8a) 

bp - Q'k = p 2 + aa - Rt - kt' , (2.8b) 

per- Qk = (p + p)a- (t + t')k + ^ , (2.8c) 

bp'-Sr' = p'p + cra' - r'f' - kk' - * 2 , (2.8d) 

together with their primed and starred versions. The Bianchi equations are given 
by 

(b -4p)*i - (8' - r0*o = -3k* 2 , (2.9a) 

(b -3p)# 2 - (8' - 2r')*i = cr'* - 2k* 3 , (2.9b) 

together with their primed and starred versions, and the Maxwell equations are 

(\>-2p)fa-(d' -r')fa = -Kfa, (2.10) 

with its primed and starred versions. Further, the GHP operators acting on 
weighted quantities satisfy the commutator relations 

[b, b'] 7] = [(f - t') 5+(t - f')3' - p(kk' - tt' + a ) - g(SR' - ff' + tf a )] 77, 

(2.11a) 

[b, 5]r? = [ - ?' b -^b' + po +ct3' - p(p' K - rV + *i) - g(cr'« - pf')} 77, (2.11b) 
together with their primed and starred versions. 

2.2. Petrov type D spacetimes. In a spacetime of type D we can fix a null tetrad 
up to rescalings (and a trivial rearrangement) by aligning the real null vectors with 
the principal null directions. Such a tetrad is called a principal tetrad. In a vacuum 
type D spacetime, working in a principal tetrad, ^0 = *i = ^3 = ^4 = follows 
and due to the Goldberg-Sachs theorem k = k' = a = a' = 0. The only non- 
vanishing GHP quantities are 

* 2 ,p, t, p'y, 

and the Bianchi identities (|2.9|) simplify to 

p*2=3p* 2 , o* 2 =3r* 2 , (2.12) 

together with their primed versions. See |16| for further identities valid for the 
GHP quantities valid in vacuum type D spacetimes. We record for later use a 



10 



S. AKSTEINER AND L. ANDERSSON 



commutation relation for {p, 0} quantities r\ on vacuum type D backgrounds. The 
following identity (and its prime) is a consequence of (|2.11b|) . (]2.8a|) and (|2.8ap *. 

[ b — ap, 3 -ar] 77 = p(5 — ar)rj — f'(b —ap)rj, (2-13) 

see |21l equation (2.5)] for a more general relation involving g- weight. 

2.3. Perturbation theory and gauge transformations. Here we give a short 
overview of gauge transformations in perturbation theory. See [471 [7] for more 
details. Perturbations of a spacetime can be understood in terms of curves in the 
space of solutions of Einstein field equations, originating at a given background 
spacetime. Linear perturbations are tangents to such curves at the origin. We 
denote the perturbation parameter by e. 

The identification of points of background and perturbed spacetime is called 
identification gauge. Introducing coordinates x a in the background, an infinitesi- 
mal transformation of the form x a — > x a + e£ a can be interpreted as changing the 
identification of points between background and perturbed spacetime. Quantities 
which do not change under these transformations are called identification or co- 
ordinate gauge invariant. A quantity is coordinate gauge invariant if and only if 
it vanishes or is a constant scalar or a constant linear combination of products of 
Kronecker deltas in the background, see |461 p. 24]. 

The Weyl scalars ^>iB transform as scalars under coordinate transformations, 

*iB^*is-e^*M + 0(e 2 ). 

In a type D background, only \P 2 is non-zero, and hence the spin scalar is the 
only one of the ^is, i = 0, . . . ,4 which fails to be coordinate gauge invariant. 

As mentioned in section |2~T1 we choose the background tetrad to be fixed up to a 
two-dimensional subgroup of the Lorentz group corresponding to boost rotations of 
the future pointing null vectors l a ,n a , and spin rotations of m a ,fh a . However, the 
perturbed tetrad has the full transformation freedom under infinitesimal elements 
of the Lorentz group 

1% n B n B + e (™ + am -), m a B -> m B + eal a , 

+ e(bm a + bm a ), n a B -> n%, m a B -> m% + ebn a , 

+ eAl a , n a B ^n%- eAn a , m% -> m% - \eQm a , 

where a, b are complex and A, are real functions (see e.g. the linearized versions 
of (|2.2p for the third line. A quantity which is invariant under these transformations 
will be called tetrad gauge invariant. For the first subset of infinitesimal Lorentz 
transformations we have for example 

*jB ->#ju + eja* i -iA, j = 0,...4;*_ia = 0. (2.14) 

A complete table for all GHP quantities can be found in [7J §5.10]. It can be verified 
that ^obt^2B and are tetrad gauge invariant. 

In the following, in order to avoid clutter in the notation, we will drop the index 
A for background quantity, unless it is not clear from the context whether a certain 
quantity is evaluated on the background. 

3. Equations for linearized gravity and electromagnetism 

In this section, we derive equations for linear perturbations of the Weyl com- 
ponents ^o, ...,-$> 4 on vacuum type D backgrounds, as well as for the linearized 
Maxwell scalars 4>o , 4>i , 4>2 on charged type D backgrounds. The gauge invariant 
fields ^>ob, ^4B satisfy the Teukolsky system TME S , for s = 2,-2, while the equa- 
tions for the tetrad gauge dependent scalars ^ib, ^33 correspond to the TME S for 
s = 1,-1, but with a non-trivial right hand side involving a gauge source function, 
cf. section 2J For the spin weight zero linearized Weyl scalar ^>2B we find a new 
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wave equation with a non-trivial right hand side. If the right hand side vanishes, 
this equation is a direct generalization of the Regge- Wheeler wave equation to type 
D. In section H] below we consider the structure of these right hand sides in more 
detail. 

We point out that while the equations for the linearized Weyl scalars as well as 
the linearized Maxwell scalars decouple in the sense that each individual equation 
involves only one of these scalars, the equations for the non-extreme spins are 
coupled via linearized spin coefficients, unless further gauge conditions are imposed. 

3.1. Weighted wave operators. As is the case for any vector bundle over (M,g) 
with covariant derivative, there is a natural generalized wave operator acting on 
sections of the bundles £{ p ' q \ The Weyl and Maxwell scalars are properly weighted 
quantities of type {2s, 0} for integer spin weights s. Since we shall be interested 
in operators acting on the Weyl and Maxwell scalars, we restrict our attention 
to the operator D p = 9 a a acting on quantities of type {p, 0}. Expanding this 
using (|2.6p and (|2.7p gives after some calculations using the commutation relations 
(l2Tia . (plTbl 



D p = 2 



Q> ~pW - p') - (3 -f'W - r') + ao> - k' k - * 2 (3.1) 

(3.2) 



+ |(/c/c' - tt 1 + pp' - a a' + 2* 2 ) 



Let oj a be the connection form in O a , cf. section [2.11 and let B a be a properly 
weighted form of type {0,0}. Then ui a — B a is again a connection form on the 
weie ted bundles and 6 a +pB a + qB a is again a weighted covariant derivative 

on the bundles £^ p ' q \ In particular, let 

B a = -(pn a - Tfh a ). (3.3) 

Modifying the covariant derivative with B a gives the weighted wave operator on 
type {p, 0} quantities 

H p = (Q a +pB a )(e a +pB a ). 

Note that Ho = H- We can now write the vacuum Teukolsky master equation for 
a spin weight s field tp^ in the form 

[H 2s - 4s 2 * 2 ]V' (s) = 0, 

cf. Bini et al. [5j §4]. 

A calculation shows that acting on a quantity of type {p, 0} we have 

H p = D p + 2pB a Q a +p{e a B a )+p 2 B a B a 

= 2(b -pp - p)(b' - pf) - 2(5 -pr - f')(5' - r') 
+ (p-2)[«/s / -o-<7 / ] + (3p-2)*2. (3.4) 
Restricting to a type D background we have 

H p = 2(b-pp - p)(b' - p') - 2(9 -pr - f')(8' - r') + (3p - 2)* 2 . (3.5) 

since k, a vanish there. 

Recalling the discussion in section \2. 11 a transforms properly under the prime 
operation. In particular, we have (Q a vY = ©a??' and hence also (<d a @ a v)' = © a ©a?7' 
for any properly weighted quantity. However, modified connection Q a + pB a does 
not satisfy this rule since B' a 7^ — B a . Instead, the operator H p has the following 
transformation rule involving rescalings. 
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Lemma 3.1. Let rj be a properly weigthed quantity of type {p, 0}. The generalized 
wave operator B p on a vacuum type D background transforms under prime as 

(B p r,y = ^'/ s B p ,(^/ 3 r)'). (3.6) 

Remark 3.2. It should be noted that thep in p denotes the weight of the quantity 
on which it acts. Hence the p' in the right hand side of (|3.6p is p' = —p, since the 
type oft)' is {-p,0}. 

Proof. Using the commutation relations and field equations on type D, one gets the 
identity 

(b' + pp> - p')Q?-p) - (b'+pr' - f)(8-r) 

= (b -p)Q?' + (p- l)p') - (3 -f')(3' + (p - l)r') + 3p* 2 . 

Rescaling the RHS by $j , using Bianchi identities (|2.12p gives the result. □ 

As we shall see, using this transformation property for E3 p , half of the equations 
for linearized gravity discussed below follow without calculation. 

3.2. Perturbation calculations. We now derive the equations for the linearized 
Weyl scalars in terms of the weighted wave operators. 

Theorem 3.3. On a vacuum type D background we have 

[Q 4 - 16*2] *os = 0, (3.8) 

[S 2 - 4* 2 ] (*^ 1/3 * 1B ) = -6*^ 3 [(J>' + 2p' - p')k b - {&' + 2r' - f)a B + 2*is], 

(3.9) 

[So + 8*2] (*2: 2/3 *2i3) = -3D B * 2 /3 , (3.10) 
[Q_ 2 - 4* 2 ] (*a X *3b) = -Q[{p + 2p- pWb - (5 + 2r - f')a' B + 2* 3 s] 

(3.11) 

[H_ 4 - 16* 2 ] (*^ 4/3 ^4s) =0, (3.12) 

Proof. First we consider the equations for linearized the Weyl scalars ^qb,^4B 
with extreme spin weights s = 2, — 2. The linearized Bianchi identities (|2.9a[) * and 
(|23a|) read 

(b' - f/)*oB = (9 -4t)# ib + 3a B * 2 , (3.13) 
(5' - t')^ob = (b -4p)*ib + 3k b * 2 . (3.14) 
Combining these identities as (b — 4p — ,o) ([3.13p — (3 — 4r — f M p.Mp gives 

[(b -ip - p)Q?' - p') - (3 -At - f')(3' - /)] *ob = 

[(b -4p - p)(3 -4t) - (3 -At - f')(b -4p)] *ib (3.15) 
+ 3[{\>-4p- p)a B - (3-4t -t')k b ]*2. 

The term involving >J>ib on the RHS vanishes due to (|2.13|) with a = 4. The 
perturbed Ricci identity (|2.8c[) reads (b — p — /o)<tb — (3 — r — t')kb = ^os and from 
(|2.12|) it follows, that also the $2 term on the RHS reduces to 3^03^2- Recalling 
the form of E3 p gives equation (|3.8p for ^os- Equation (|3.12p for follows from 
this after applying a prime and using (|3.6[) . 

As we shall see, the corresponding wave equations governing the linearized Weyl 
scalars ^ib,^3B with spin weights s = 1,-1 do not decouple in the sense that 
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other perturbed quantities than the given linearized Weyl scalar are involved. The 
linearized Bianchi identities (|2.9b[) * and (j2.9bp ' read 



(3.16) 



(b'-2p')*i fl = |(9-3r)* 
-(3'-2t')*ib = -|(b-3p)f 2 J , 
Multiplying both equations by \E , 2 ~ 1/ ' 3 anc ^ usm g the Leibniz rule gives 

(b' - P ')(^ 1/3 *ib) = ^|(3-2r)*f|^, (3.17) 

-(3' - r'){^- 1/ ^ 1B ) = -^|(b-2p)^ /3 |^. (3.18) 

Here we have made use of the fact that ^ 2 1 can be moved inside the {}b brackets 
in view of the background Bianchi identities (|2.12|) . 

Combining the above identities as (b ~2p~ p) (|3.17p + (3 — 2r — -f ') (J3H5J) , we get 



(b -2p - p)(b' - p') - (8-2t - f')(3' - t<) 



| [(b -2p - p)(3 -2r) - (3 -2t - f')(b -2p)] ^ /3 | 



which gives 

[H 2 - 4* 2 ] {^ 1/3 <S> 1B ) = 3| [Q> -2p - p)(3 -2r) - (3 -2r - f')(b -2p)] *\ fz J , 

(3.19) 

using (|3.5p . Expanding the right hand side of equation (|3.19p , leaving off a factor 
of 3, we have 

[Q> -2p - p)(3 -2r) - (3 -2t - f')(b -2p)] B ^/ 3 + 

+ [(b -2p - p)(3 -2r) - (3 -2r - f')(b -2p)] ^ /3 ] fl . 

The second term vanishes due to the commutation relation (|2.13p with a = 2, but 
for the first term, we must use identites valid off the background. Since ^ is a 
{0, 0} quantity, we find 

[(b -2p - p)(3 -2r) - (3 -2r - f')(b -2p)]a*& 3 = 
= [-f'b-Kb' + p3+er3'-2Tb-2(bT) + 

+ 2p3+2(3p) - (2p + p)(3-2r) + (2r + f')(b -2p)] B ^ (3.20) 

= [-2/tp' + 2ot' - 2(b r) + 2(3 p) + 2rp - 2f 'p^*^ 3 

= -2^ 3 [(b' + 2p' - p')K B - (3' + 2r' - f)a B + 2*1*], 

where we used the commutation relation (|2.11bp in the first step, background 
Bianchi identities in the second step and the Ricci identities (|2.8ap ' and (|2.8aD 
in the last step. Applying a prime and making use of (|3.6p gives (|3.1ip . 

Finally we consider the spin weight linearized Weyl scalar ^2B- Recall that 
*2 is non- vanishing in a type D spacetime. For this reason, it is convenient in the 
calculations to leave some expressions as {} B brackets. As in the previous cases, 
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we start with linearized Bianchi identities. From (|2.9b[) ' and t|2.9b[) * / . we get 

|(t>' - 3p')* 2 | = (3-2r)* 



|(8'-3r')* 2 J = -(b-2p)* 3B . 



(3.21) 



— 2/3 

Multiplying both equations by * 2 ' and using the Leibniz rule gives rescaled 
equations 



3(b' - P ')*l /3 
| 3 (9'-r')^ /3 



v&- 2/3 (3-2r)vI/3 i j, 



= -* 2 2/3 (p-2p)*3B. 



(3.22) 
(3.23) 



Here we used the fact that * 2 2 ^ 3 can be moved inside the {} B brackets because of 
the background Bianchi identities (|2.12[) . 

To find a wave equation for *2B, we consider the combination fb — p) (|3.22p + 
(5— fQ lEggg)) , which gives 



(b-p)^3(b'-p')* 2 /3 



f) 3(9'-r')< 3 



(b -p) (*~ 2/3 (3 -2t)* 3B ) - (3 -f') (*" 2/3 (b -2p)* 3B ) . (3.24) 



Using the identities (b — p)(* 



-2/3,^ _ lT ,-2/3 



*2 Q> - 2 / - P)^ and ( 8 -^')(* 



-2/3, 



*2 2/3 ( 8 - 2r ^ f which follow from (|2~T2")) . the RHS vanishes due to ([2TT5f with 
a = 2. The operators on the LHS can be put into the {}b brackets because of the 
background Bianchi identities (|2.12[) . This gives the identity 



□o + 2* 2 ]* 



1/3 



= 0. 



Expanding the {}b bracket in equation (|3.25p gives 



= 3j[Q + 2* 2 ]* 



1/3 

2 



□ + 2* 2 



*2 ' *2I 



-2/3 
-2/3 



□ B + 2* 2B 



* 



1/3 



*2 ' *2B +3D B * 



,1/3 



and hence p.!0[) . This completes the proof. 



(3.25) 

(3.26a) 
(3.26b) 

(3.26c) 

□ 



For future reference, we state the following equations which were used in the 
proof of theorem [ 

Corollary 3.4. 

[S 2 - 4* 2 ] (* 2 T 1/3 * 1B ) =3U(P -2p - p)(d-2r) - (<5-2r - f')(P -2p)] * 2/3 J 



[□0 + 2*2]* 



2 /3 



(3.27a) 
(3.27b) 



[□2 - 4*2] (* 2 r 1 *3B) = 3*~ 2/3 j [(P-2p - p)(d-2r) - (<5-2r - f')(J>-2p)]* 

(3.27c) 



2/3 
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The rescaled Bianchi identity p,18[) has the same form as the perturbed Maxwell 
equation (|2.10[) on a charged type D background {<f>\A 7^ 0) 

j(b-2p)0ij =(& -t'Wqb. (3.28) 

Therefore the decoupled electromagnetic perturbation equations follow immedi- 
ately, and we have the following result. 

Corollary 3.5. In a charged type D background the Maxwell components (j)iB,i = 
0, 1, 2 fulfill the equations 

[0 2 - 4* 2 ]0 OB = 2 j [(P-2p - p)(d-2r) - (<J-2t - f')( J &-2 /9 )]0 1 | (3.29) 
{ [{P-P - P)(P' - V) - (0-r - f'){d' - 2r')] 0i } = (3.30) 

[H' 2 - 4* 2 ]0 2B = 2| [(J>-2p - p){d-2r) - (5-2r - f')(J>-2p)] Vij (3.31) 

For a charged type D spacetime, the background Bianchi identities (|2.9p include 
a term involving 0i0i and hence the simple rescaling used above for vacuum type D 
backgrounds does not apply. Instead one can use the background Maxwell equations 
(|2~TTJ|) . namely 

p0i = 2/90i, o0i = 2t0i, 
as in [18]. But as the 5' 2 rescaling is singular for the limit of flat background, the 
01 rescaling does not work for the limit of uncharged background. 

In the special case of a test Maxwell field on an uncharged background, rescaling 
by ^ 2 becomes possible and the equation for 0i B reduces to the Fackerell-Ipser 
equation [TS] 

[□ + 2* 2 ](*- 1/3 0i B ) =0, 
while the equations for 0os and 2B become the spin s = ±1 TME S . 

4. Gauge source functions 

In this section we consider the equations for the gauge dependent quantities 
^ib, ^2B, ^3B in more detail. 

4.1. Gauge source functions for the Einstein equations. In [23], Friedrich 
derived a frame based, symmetric hyperbolic system for the Einstein- Yang-Mills 
system. We specialize to the vacuum case and set the conformal factor = 1. 
Then, the result of [23] gives a symmetric hyperbolic system for a set of unknowns 
consisting of a null tetrad (or spin frame), the spin coefficients and the Weyl spinor. 

Starting from a system of equations involving the tetrad and the curvature com- 
ponents as variables, he identified the gauge source functions for this system. These 
are, letting (x^) be coordinates on M, 

and 

F ab = V a (y a e b a )e c b9bci 

where (e a ) is a null tetrad, cf. [23J equations (2.6), (2.13)], see also [24"] . 

For the frame based hyperbolic system considered by Friedrich one may freely 
specify the gauge source functions as functions of the spacetime coordinate, the 
tetrad, the connection coefficinets, and the Weyl tensor components, i.e. 

F» = F»(x a ,e«,T a bl ,W abc0 ), (4.1) 

F ab =F ab (x a ,e«,Tl,W ablt) ) 7 (4.2) 
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without changing the principal part of the resulting symmetric hyperbolic system, 
see the discussion in [23], in particular [23] p. 1462]. As we are in a geometric 
situation where it is natural to adapt to a specific background geometry and to use 
a GHP weighted tetrad, it is convenient to consider the following modified gauge 
source functions. They differ from the expressions given by Friedrich by lower order 
terms, which do not change the principal part of the resulting reduced system. 

To define the coordinate gauge source function, fix a background metric g a b on 
M, with Levi-Civita derivative V a and let V a be the tension field for the identity 
map (M,g a b) -> (M,g ab ) defined by 

V a ta = <7 cd (V c - V C )U 

holds for any 1-form £ a , see [2J for details. Then a gauge source function for the 
coordinate degrees of freedom can be given by the equation 

F>* = V^. (4.3) 

Further, let O a be the weighted GHP covariant derivative and let (e„) be a weighted 
GHP tetrad. A gauge source function for the tetrad degrees of freedom can be given 
by the equation 

F ab = Q a (Q a e b a )e c b9bc . (4.4) 
These are the expressions which we shall consider below. 

4.1.1. Gauge source functions for the Linearized Einstein equations. Consider lin- 
earized perturbations around vacuum spacetime (M,g a b). The work of Friedrich 
on hyperbolic reductions carries over immediately to the linearized vacuum field 
equations. Thus we may consider linearized frame based systems with unknowns 
consisting of the linearized tetrad, the linearized spin coefficients and the linearized 
Weyl scalars. The reduced system is extracted by specifying linearized coordinate 
and frame gauge source functions F B , F Ba b which may be specified freely as func- 
tions of the unknowns which are linear as functions of e Ba , F Bbc , WsabcB- Thus, 
with this restriction, we may consider gauge source functions 

F£=F£(x a ,e% a ,T% bc ,W B abt*), 

Fsab — F Ba b{x a , e Bal ^Bbti ^Babco)- 

4.2. Linearized Weyl scalars and gauge. In the following discussion it is in 
some steps convenient to use a 5 to denote first order linearized fields. In particular, 
5g a b = h a b and denote the resulting perturbations in geometric fields defined in 
terms of g a b by e.g. 8R a b- We have 

SRab = — -^h a b — Ra C b d h c d + V^Uft), 

where, letting h = g ab h a b, and working in a coordinate system (x^), 
v a = Vph^ - ^V a h = g Pu 5T%. 

The vector field v a defined by this expression is precisely the linearization of the 
tension field V a around g a b (playing the role of the background metric g a b above) . 
Thus v a is the appropriate coordinate gauge source function for linearized pertur- 
bations h a /3 of g a /3 and the gauge condition which corresponds to (|4.3[) is given 
by 

/ = Fg. 

The standard harmonic gauge, also known as deDonder gauge, with respect to the 
background metric is given by the condition F B =0. In this gauge, the linearized 
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Einstein equations in terms of the linearized metric h a {, take the form of a wave 
equation 

Oh a b + 2R a c b d h c d = 0, 
where □ = V C V C is the covariant d'Alembertian. 

Let (M, g) be a vacuum type D spacetime. Working in a principal null tetrad, 

1 /3 

let ip = ^ 2 . Recall that ^ is of type {0, 0} and is thus a well-defined function on 
spacetime. In notation used in this section, the condition that the right hand side 
of (|3.10p vanishes, i.e., 

□ B * 2 /3 = 0, 

takes the form 

(6n)ip = 0. 

We have 

(6D)il> = -h ab V a V b iP - v a V a il>. (4.5) 
We can view this equation as specifying part of the coordinate gauge degrees of 
freedom. In the Schwarzschild case, working in a principal tetrad, ijj is real, and 
hence (|4.5[) specifies one component of v a . In the general case, -0 is complex, while 
v a is real. Taking the real and imaginary parts of (|4. 5|) gives two real equations for 
v a . 

In order to analyze this equation in the Kerr case, it is convenient to calculate in 
a coordinate system and tetrad which is non-singular on the horizon. A tetrad in 
the ingoing Kerr coordinate system (also known as ingoing Eddington-Finkelstein 
coordinates) was described by Teukolsky §5]. In this tetrad, the components 
VaV 1 are non- vanishing on the horizon. 

As mentioned above, it is compatible with the well-posedness of the reduced field 
equations to allow the gauge source functions to depend on the Weyl scalars. Thus 
we may also consider gauge conditions of the form 

- 3Db>F 2 /3 - 6<F 2 /3 <F 2 b = 0, (4.6) 

which leads to the wave equation 

(□ + 2* 2 )(* 2 r 2/3 * 2B ) =0, (4.7) 

for ^2B- Thus, in the gauge given by ()4.6[) . >F 2 2 ^ 3 ^2B satisfies the Fackerell-Ipser 
equation. This substantiates the discussion in the work of Crossman and Fackerell 
[IH 118) . In the Kerr case, calculation shows that the gauge source function given 
by (|4.6|) will have terms depending on 1/a, and hence this gauge condition behaves 
in a singular manner in the Schwarzschild limit. Equation (|4.7p is not known to be 
separable or admit a symmetry operator, see however [181 p. 617]. This discussion 
shows that in the rotating case, also a generalized harmonic gauge condition leading 
to a homogenous wave equation 

□(*a 2/3 *2 B ) = o, 

(which admits symmetry operators) is compatible with a well-posed reduced system. 
Next we consider the phantom gauge condition. A calculation shows that 

e a m b QJ b = 2[ P 'k - Q'ci - p'k + ?a + 

Restricting to a type D background this expression vanishes, and the first order 
linearization gives 

F B i m ■= ije°m 6 e o Z 6 j =[0>' -p')n B -{b' -t>b+*ib]. (4.8) 
Now we can write the equation for in the form 

[Q 2 - 4* 2 ] (^^ib) = -^I'aWbuu + 2 P 'k b + 2r'a B + *ib]. 
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Thus, Chandrasekhar's phantom gauge condition can be written in terms of the 
gauge source function as 

F B im = -2p'« B - It' o B - tf 1B . 

In view of the discussion above, this form of the gauge source function is compatible 
with a hyperbolic system for the linearized Einstein equations. The equation for 
can be handled along the same lines. 

5. Gauge invariant equations on Schwarzschild background 

Price [55] has shown that for linearized gravity on a Schwarzschild background, 
r 3 3m^2B describes odd parity perturbations. He used the special coordinates of 
Newman and Penrose |36l p. 572] for the perturbed spacetime, which can be un- 
derstood as a coordinate gauge. Price expressed the perturbed spin coefficients, 
which occur in the ^2B equation, in terms of perturbed metric coefficients. These 
coefficients are real in the odd parity case and therefore cancel. Price then used the 
definition of CJm v E , 2B hi terms of perturbed Riemann tensor components, to relate 
it to the perturbed metric coefficients. He showed that these coincide up to a time 
derivative, in Regge- Wheeler (RW) gauge, with the RW variable Q for odd parity 
perturbations. 

Starting from (l3.2Gcl) . 

* 2 ~ 2/3 *2b) = -3M 2 /3 , (5.1) 



□ + 8* 2 

we rederive the known result, that the imaginary part satisfies the gauge invariant 
RW equation [35] , which reduces to the result of Price in RW gauge. We also show 
that the real part gives the gauge invariant Zerilli equation. Price was not able to 
derive this equation, possibly due to his choice of special coordinates. 

The conventions of Martel and Poisson [M| will be used in this section. Indices 
a,b, . . . for coordinates t, r and A, B, . , . for coordinates 9, cp. This notation differs 
from the convention used by Regge- Wheeler and Zerilli in that A'mp = -Krw - 
1)G. Metric perturbations are denoted g^ = gAfiv+Pfiv, since is used for some 
spherical harmonic decomposed components. It should be noted that we still use 

the signature (H ) while Martel and Poisson use the signature ( — h++). Using 

Schwarzschild coordinates g a b = diag(/, — — r 2 , — r 2 sin 2 #), where / = 1 — ^jf- 
and ^2 = —AI/r 3 , the Kinnersley frame reads 

l a = (/-\ 1,0,0), n Q = i(l,-/,0,0), m a = -L- ( , 0, 1, ^-). 
v ' 2 V ' V2r sintr 

With V = *2 /3 > the RHS 01 (EH> reduces to 

= ( - iTdn - (V) - p r T£, + g rr (d r p,n)d r ^ 

Further simplifications will occur in odd and even part, which are investigated in 
the next sections. 

To relate ^>2B to Regge- Wheeler and Zerilli variables we do a calculation analo- 
gous to that of Price [421 Appendix D]. Starting from 

-2^ 2B =R AaM s{l a n l5 Vn S ) B ~R AaMS {l a n^m^m S ) B 

3 8 8 i-8 ( 5 ' 3 ) 

+ RBaBjs{l a n p Pn )a~ RBaB 1 8(l a n p rn l m ) A , 

we get 

-2* 2B = RBtrtr + RBrt9<p + —Pafil"^ . (5.4) 
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The second term is purely imaginary and known from the calculations of Price. 
The perturbed Riemann tensor is related to metric perturbations via 



3~f;a5 + Pa6;/3~f Pa*y;@5 Pft6;a~f 



With these relations one can check explicitly that 3m x f , 2B corresponds to odd parity 
and d\^2B corresponds to even parity perturbations. 



5.1. Imaginary part and Regge Wheeler equation. The odd parity metric 
perturbation can be expressed as [3H eq. 5.1 - 3] 



Im 
AB ' 



where X l B n and X l ^ B are vector and tensor spherical harmonics. It follows that 
trp M! y = and after a short calculation (|5.2p reduces to 

Db^j = 0. 

3m*2S corresponds to odd parity perturbations and 3m (|5.ip reduces in this case 
to 

[□ + 8* 2 ](* 2 ; 2/3 3m* 2i j) = 0. 

Introducing tortoise coordinates r* by <9„ = fd r cancels 2^2 from the potential, 
a rescaling by r cancels first order d r terms and we are left with a gauge invariant 
Regge Wheeler equation 



1(1 + 1) 



-f 



6M 



(r 3 3m*2B) = 0. 



To relate this to the Regge Wheeler variable Q, we look at (|5.4p . For odd parity 
perturbations it gives {R° B dd )trtr = = p°A d l a n p and with ([53]) . 



r 2 sin 6 



1(1 + 1) 



(for convenience we suppress spherical harmonics and the related indices). It follows 
that just the imaginary part contributes to the perturbations. We now have 



1(1 + 2)! 



4r(7-2)! 



— ho t ~~ — ho 

2 ' r 



j .i 



(* + 2)! 
4(Z - 2)! 



-iodd 



which is the gauge invariant variable of Moncrief [33] and in RW gauge reduces to 
the result of Price [12] • 



5.2. Real part and Zerilli equation. The even parity metric perturbations can 
be written [Ml eq. 4.1-3] 

Pab=2_^ h ab Y . PaB=2^Ja Y B J PAB = r ^ K UabY + G Y ABi 
Im Im Im 

(5.6) 

where Y lm ,Y B m and Yj^ B denotes the even parity scalar, vector and tensor spherical 
harmonics. U\c^2B is not gauge invariant, but transforms as 



2B 



3Ai\ r 

- -j-r- 
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For some coordinate gauge transformation 



£ M . Therefore we rewrite the 



wave equation for the gauge invariant quantity, using Appendix [Al 



3M 



2B 



r d 9te* 2 B + — [K + -G 



X 



(5.7) 



where A = 1(1 + 1). The ^2B equation f|5 . 1 1) for even parity perturbations reads 
A „6M" 



d 2 t-dl+f~ 2 -f , 



3Mfr\3 B r 



r 
3M 



*2B = 



A 



at - % m + /-j - / 



.6M 



A" 



A 



G* 



The perturbed wave operator term (|5.2p does not vanish, but gives 



— Jr 



(ZMf f 2 



M 



-A 



(5.8) 



/a r A . 

(5.9) 



Expanding (|5.4[) for even parity perturbations gives after some calculations Rsrt 
and 

JKetf 2 B =M(f- 1 hu ~ fhrr) 



3 r 



+ <9 t 2 /i r , 



2m 



<9r(/ ^ti - //irr) J^d t h tr ~ 2d't r h tr 

f J 

(5.10) 

The whole equation is now given in terms of metric perturbations. To compare 
the result to others, we express all metric components in terms of gauge invariants 
(denote with a tilde according to the conventions of Martel and Poisson). Using 
the equations of Appendix [21 ^2B takes the form 



2B 



Ar 

T 



K 



2/ 
A 



fh rr — rd r K 



(5.11) 



where A = 1(1 + 1), \i = (I - l)(l + 2) = A - 2, A = \i + 6M/r. The RHS of lETSll 
can also be expressed in terms of gauge invariant quantities 

RHS = -/!f5 2B + ^A-. (5.12) 

The first term cancels the RW potential on the LHS! Now we rescale (|5.8p by A -1 
(which depends on r), 



(A-^as) = 



3M/A A _ X 



2r 2 



A- i A-2(9 r ,A- i )(a r ,* 2B ) 



(5.13) 

A straight forward but tedious calculation shows that the RHS can be written as 

(5.14) 



^A^A - 2(d r ,A^)(d r ^ 2B ) = f^A~^ 2B , 
2H r J A 



the new potential term on the LHS simplifies to 

12M/ / 6M 2 



HdlA- 1 ) 
and we finally end up with 



A 2 7 



3Mfi — fir 



1J2 



(5.15) 

= 0. 

(5.16) 

This is the gauge invariant Zerilli equation. The relation to Moncrief's gauge in- 
variant variable is simply Q cvcn = 4A~ 1, I'2B- 



12M/ f6M 2 , u 2 r , M -r 

3/xM + — 



A 2 r 4 V r 



6M 



£ + 1 



(A" 1 * 



2Bj 
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Appendix A. Even partity perturbations in Schwarzschild coordinates 
We used equations of the unpublished appendix of [34 . which is available as 



arXiv:gr-qc/0502028. For convenience we repeat the required results of appendix 



C. 

The even-parity metric perturbations (|5.6p transform under even parity gauge 

£ ^Im^lm ^Im^hn ^my?mj 

9 . 2Mf^ PI 9, 9, 2AT 

3, 2M, 9, . 

Sh rr = -2frSr ~ —^r, Sjt = - ft, 

9 2 2f A 

#r = -^-f-fr + -f, = Kr + "j f, 

or r r r z 

Martel and Poisson extracted the following gauge invariant quantities 

- , 9 9 2M , d 2 „ r - 3M 9 ^ 
h tr = h tr - - Jt - - Jr + ^j Jt + r —G + J Qt ' 

h -h o 9 ■ 2M ■ 2 9 2 r , 2r - 3M 9 
ar r z / ar^ j or 

r or 2 

and with these, the vacuum field equations are 

0- ° 2 K 3? " - 5M 9 K \- ——h (A + 2)r + 4M - ^ ^ 

9r 2 r 2 / 9?* r 9r rr 2r 3 rr 2r 2 / ' 

0— —K— — — K-^ — h —h 

dtdr r 2 f dt r dt " 2r 2 * r ' 

n 9 2 T y (r - M)f 9 - 2/ 9 - /9- Ar + 4M - / 2 - ,i/ - 

at r z or r at r or if 5 r z 2r 

n _9~ 9 - 19- 2M- 

9t 9r j 9r rv 

9 - 9 - . 9 - r-M- (r-M)f- 

at or or r z j r A 



rr • 



9 2 - 9 2 - 9 2 - 1 9 2 - 9 2 - 

= _ — hrr + 2—h tr - — h tt - j—K + f—K 

2(r - M) 9 - r - 3M 9 - (r - M)/ 9 - 2(r - M) 9 - 

H 77 'KL' l tr 77 — -^-h tt ~ -^-h rr H = TJ - -^ 

r z j at r z j Or r l or r z or 

Xr 2 -2(2 + X)Mr + 4:M 2 ~ Ar 2 - 2fxMr - AM 2 - 

-"■tt ^ — i h r 



2r 4 f 2 2r 4 



0=jh u - fh r 
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